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Abstract Coronal mass ejections (CMEs), which dynamically connect the solar 
surface to the far reaches of interplanetary space, represent a major manifestation 

^j ■ of solar activity. They are not only of principal interest but also play a pivotal 

«, \ role in the context of space weather predictions. The steady improvement of 

both numerical methods and computational resources during recent years has 
allowed for the creation of increasingly realistic models of interplanetary CMEs 
(ICMEs), which can now be compared to high-quality observational data from 
various space-bound missions. This review discusses existing models of CMEs, 
characterizing them by scientific aim and scope, CME initiation method, and 
physical effects included, thereby stressing the importance of fully 3-D ('47r') 

Ql{ ■ spatial coverage. 



Keywords: Coronal mass ejections: initiation and propagation, Coronal heat- 
ing theory, Magnetohydrodynamics, Solar Wind: disturbances. 



o. 

£3 ■ 1. Introduction 

CO 



Coronal mass ejections (CMEs) are one of the most spectacular manifestations 
of solar activity. Their mass and energy output, insignificant as they may seem 
when compared to the Sun as a whole, are still able to severely distort plane- 
tary magnetospheres even after they have been diluted considerably during an 

qq ■ expansion phase covering a distance of one AU or more, thus bringing the solar 

influence to the outermost reaches of interplanetary space. 

^sO ■ CMEs also relate to many othe r fields of solar and stellar physics . They are 



^C) \ intimately linked to flares (e.g. IShanmugaraiu. Moon, and VrsnaH 1201 lL and 

references therein), and often give rise to strong shock fr onts, which act as 
efficient accelerator engines for solar energetic particles (e.g. iReamed . Il999l ). It 



has also been recognized that CMEs play an important role for the restruc- 
turing of the Sun's global magnetic fi eld over the course of the solar cycle 
(jSchwadron. Owens, and Crookerl . 120081) . 



Apart from these motivations to understand the CME phenomenon from a 

C^ ■ purely scientific perspective, it has also become apparent in recent years that 

the increasing complexity of modern technology, in particular communication 
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infrastructure, has also led to an increased vulnerability of this technology to the 
adverse effects of space weather. The resulting commercially-driven demand for 
high-precisi on forecasting tools has given a continuing boost to research efforts 
in the field ( Bakeri l2002t IPulkkinenl 120071 and references therein) . Besides the 



ensuing increased public interest in the subject, timely advances in computing 
hardware and numerical algorithms have allowed for the development of sophisti- 
cated large-scale models of CME dynamics, which are the subject of this review. 
The paper is structured as follows. After this introduction, Section [2] will discuss 
ways in which the vast spectrum of existing CME models can be categorized, 
with special focus on their specific strengths and possible weaknesses. Section [3] 
describes how different types of models can be validated against obscvations, and 
Section 0] presents a concluding summary and a modest suggestion for further 
action to improve on the comparability of results from different models. 



1.1. What makes CME modeling a demanding task? 

To realistically model the evolution of a CME is a very demanding task for 
several reasons: 

1. The CME phenomenon spans vast temporal and spatial scales. Even disre- 
garding the microphysics involved in, e.g., flare reconnection (some 10~ 8 s, 
several meters) best described with kinetic approaches, a CME's life cycle 
from pre-eruption (when it extends over a small fraction of a solar radius 
Rq) to interplanetary expansion and finally merging with features like global 
interaction regions at the far reaches of the heliosphere ex t ends o ver some six 
orders of magnitude in space and time, see iForbes et al\ ( 20061) and in par- 
ticular Figure 1 therein. For this reason, existing fluid-based models usually 
specialize on either initiation/eruption, interplanetary expansion, or interac- 
tion with corotating interaction regions, planetary magnetospheres, or other 
CMEs. 

2. Even single CMEs show great variety in their morphology. Although during 
solar minimum, many of them originate from streamer blowouts and often ex- 
hibit the famous three- part "light bulb" structure consi sting of a bright, semi- 
spher ical front enclosing a dark cavity and a bright core ( Illing and Hundhausenl .[ 
Il986l) . the situation changes towards solar maximum; CME events then usu- 
ally originate from active regions and tend to exhibit a much more irregular 
structure, which can differ markedly from the three-part textbook config- 
uration. Even for prominence-related CMEs, the appearance of the three- 
part structure depends on the prominence location (jCremades and Bothmerl . 
200J). The fraction of CMEs that canno t be classified into the subgroup 



of 'magnetic c louds' (jBurlaga et all . 119811 ) was estimated to be near 2/3 by 
Gosling ( 19901) . albeit with a marked variation throughout the solar cycle 
from almost no ne d uring some years of minimum up to ~60% at maximum 



( Li et q/.l . l201ll ). See lKilpua et al\ (J2011I) for a recent review on magnetic cloud 
models and related multipo int observations. In their sample of almost 1,000 
CMEs. lHoward et ali (J1985I) were able to identify as many as ten morpholog- 
ical classes, and even this list is probably still far from exhaust ive. On top of 
this, the large angular width of CMEs (sa 50° on average, see ISt. Cvr et al. 
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(2000)), combined with their high rate of occurrence especially during solar 



maximum, makes interaction among them likely, thus giving rise to an even 
wider spectrum of possible morphologies. 

3. Attempts to reproduce observed events tend to be severely under-determined 
from observational side. To fully constrain a well-posed magnetohydrody- 
namic (MHD) model, initial and boundary conditions for density, pressure, 
and the velocity and magnetic field vectors must be specified; yet the initial 
(pre-eruptive) conditions are poorly known and usually rely on surface mag- 
netograms and coronagraph images to constrain the magnetic field structure, 
possibly complemented by in situ observations at individual locations. 

4. CME propagation is an inherently three-dimensional (3-D) process. Even if 
the CME initially exhibits some form of spatial symmetry, its ensuing expan- 
sion in a complex magnetic environment will in any case break that symmetry, 
rendering modeling approaches with implied spatial symmetries problematic 
(see Section ET231) . 



2. Model classifications 

Both the complexity and the enormous range of scales covered by the CME 
phenomenon make it necessary to devise separate models to address the different 
stages of CME evolution, most notably the onset and eruption, the ensuing phase 
of propagation and expansion, and a possibly following interaction with other 
structures (planetary magnetospheres, other CMEs). 

2.1. Modeling CME onset... 

Since the kinetic energy of a CME is in most cases greater than what the 
photosphere can provide during the timespan of the eruption, all current models 
for the initiation phase assume a slow accumulation of energy, followed by its 
rapid release due to loss of equilibrium. Existing analytical models for this early 
phase are few in number and have to rely on simplifying assumptions to keep 
the calculations manageable. 
In the 2-D ideal MHD model of lForbes and Isenbera (1991), the magnetic energy 



of an initially stable filament is slowly increased by converging advection of 
additional flux until a current sheet forms. Reconnection below the filament 
can then lead to a 'catastrophe', i. e., a sudden loss of mechanical equilibrium 



which causes the filament to erupt. iLin et all (|1998f ) extended this model and 



considered a finite curvature along the invariant direction (thereby transforming 
the infinitely long cylindrical filament into a force-free toroidal flux rope around 
the Sun). While both models predict eruption as soon as the stored magnetic 
energy exceeds a certain threshold, the now finite curvature force was shown to 
yield a qualitatively different behavior, which now favors eruption for large flux 
ropes. 

Early numerical approaches to the simulation of solar eruptions relied on re- 
latively simple 2-D arcade geometries, in which a forced shearing motion of the 
structure footpoints was employed to generate slow CMEs by driving the system 
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past a critical point in its parameter space (e.g. iMikic. Barnes, and Schnack , 
19881 : ISteinolfsonl . Il99lh . 

The currently considered scenarios for CME outbreak largely fall into the fol- 
lowing groups: 

1. The mass loading model relies on the notion that within a prominence, the 
equilibrium between magnetic tension and the gravity force of the mass above 
the structure allows for the slow accumulation of more mass, which merely 
cause s the prominence to bend downwards under the in creased weight of its 
load (IFong. Low, and Far! 12002c I Zhang and Lowl 120041 ) . Eruption can occur 
when some restructuring of the magnetic field causes part of the material to 
suddenly drain away, causing a loss of equilibrium. These authors distinguish 
'normal' and 'inverse' prominences based on their orientation with respect 
to the background field, and find that the normal vari ant is more l i kely t o 
erupt. This finding has been confirmed numerically bv IChane et all i 20061) . 
see Section |3~T1 

2. Flux cancellation models (also known as 'catastrophe models') start with a 
(usually twis ted) arch-shaped flux rope w hich undergoes reconnective flux 
cancellation ( Martin. Livi. and WaneLll985 ) at its neutral line. This increases 
both the tube twist and its magnetic pressure at the expense of the overlying 
field, thus moving the equilibrium position to greater heights until eruption 
occurs because no further neighboring equilibrium exists. Notable numerical 
investigations o f this process were c arried out, in incre asing order of com 



Plexity , by, 
(|2004h . 



e.g. 



Amari et all ( 2000J) , iLinker et all (J2003I ) , and iRoussev et al. 



3. Tether cutting (jMoore et aUuOOll) is very similar to flux cancellation, except 



maybe more impulsive and at slightly larger photospheric heights, cf. IChenl 
( 20111 ). It starts with a single closed bipole (essentially a magnetic arcade) 
consisting of a twisted sigmoidal core and an 'envelope' of less twisted field 
lines. If the core axis is suitably aligned with respect to the neutral line, 
reconnection below the axis will weaken the tension of the envelope. This 
causes the core part to move upwards, thereby stretching the envelope field 
and leading to even more reconnection in the current sheet thus formed below 
until confinement becomes too weak to prevent the core from erupting. The 
model agrees well with x-ray images from the Soft X-ray Telescope (SXT) on 
board Yohkoh, but is deemed implausible o n the basis of energy considerations 



-, 



(jAntiochos. DeVore. and Klimchuk, 1999). 



4. The breakout model (lAntiochos. DeVore. and Klimchukl . Il999i) starts with a 
quadrupolar arcade featuring an external X-type null point. The arcade is 
sheared and/or twisted by footpoint motions until the null is deformed into a 
current sheet at which reconnection sets in. This removes some of the overlying 
flux, which widens the current sheet further, thus leading to a runaway process 
which eventually causes the arcade to erupt. This model's crucial feature is 
that the ejected plasmoid is topologically detached from the arcade right 
from the b eginning, which is a favorab le condition to avoid the Aly-Sturrock 
constraint (JAM Il984t ISturrocH Il991 ), see below. Corresp ondin g MHD sim 



ulations have been carried out by IMacNeice et al.\ (J2004I ) and ILvnch et al. 
(l2004h 
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5. Finally, two different variants of the MHD kink instability h ave been proposed 
as pos sible initiation mechanis ms. Using the analyt i c field bv lTitov and Demoulin[ 
( 19991 ) as initial configuration, iTorok and Klieml ( 20051 ) showed that if this 
flux rope is twisted beyond a critical value, the twist is partially transformed 
into writhe, and a helical MHD ki nk instability sets in, as was already hinted 
at by an earlier stability analysis (JHood and Priestl 1 1 9 8 lh . This may or may 
not l ead to an eruption, depending on the strength of th e overlying field (e.g. 
Fanl . l2005l ; lRachmeler. DeForest. and Kankelborgl . 120091 ). The kink instability 



thus not only addresses and explains the sudden release, but also the previous 
storage of magnetic energy, and can furthermore explain some relat ed obser- 
vations such as soft x-ray s i gmoid s (|T5rok. Kliem. and Titovl . 120041 ) . 



Second, iKliem and Torokl (120061 ) analyzed the stability of a torus-shaped 



current loop against radial perturbations and the possibly ensuing torus insta- 
bility, which they described as " a lateral kink instability distributed uniformly 
over the ring." They found that this instability may trigger the onset of a 
sclf-accclcrating expansion along the major radius that can explain both the 
eruption of a flux rope (when treated as the upper half of a torus whose 
lower half is submerged below the photosphere) and the ensuing acceleration 
behavior of different types of CMEs, notably slow and fast ones. 

Common features of the last four models are the formation of a twisted flux 
rope and the generation of a current sheet below the flux rope. What sets the 
first and last model apart from the others is that they rely on purely ideal 
MHD effects, and do not involve any reconnection (at least not as a trigger 
clement). Reconnection has been viewed as a convenient means to circumvent 
the Aly-Sturrock constraint according to which any process which entirely opens 
a force- free field to infinity would in fact increase the overall magnetic energy and 
therefore could not drive an eruption. However. iForbes et al\ ( 20061) enumerate 
several other avenues to avoid the constraint (such as fields being non-force- 
free, not sim ply connected, or containing field lines d i sconn ected from the Sun). 
Additionally, iRachmeler. DeForest. and Kankelborgl ( 20091 ) used a Lagrangian 
simulation scheme with no (numerical or physical) resistivity whatsoever to 
demonstrate that reconnection is not necessary to drive fast CME eruptions. This 
theoretical finding is corroborated by the small but non-zero number of observed 
fast CMEs showing no sign of flare association ([Marque. Posner. and Kleinl . 
20061 ). More details about the initiation and early propagation phase of CMEs, 
especially of those origi natin g from helmet stream ers a nd prom i nences , can be 
found in the revie ws by [Cowl (1 19961) . IForbes! (|2000l ). and lVrsnakl (J2008[) . as well 
as in Cha pter 8 of iHowardl (l201lh. 

Recently, iLin. Gallagher, and Raftervl ( 20101 ) analyzed the relevance of some of 
the above-mentioned mechanisms for the kinematic evolution of two CME events 
and found the breakout and catastrophe models to yield the best agreement 
with their derived time profi les of he i ght, v elocity, and acceleration. This is 
reminiscent to speculations bv IHowardl ( 20111) suggesting that "...it is likely that 
different types of CMEs are best described by different models. Indeed, it is 
possible that most, if not all of the models [...] may be appropriate to describe 
some CMEs under certain conditions." 
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2.2. ...and propagation 
2.2.1. Trajectory mapping 



For the purpose of space weather forecasts, the three most crucial quantities 
to be delivered by a model are the CME's trajectory, travel time, and geo- 
cffcctivencss. To a first approximat ion, expansion i s radi al, implying that only 
halo CMEs are likely to hit Earth. ISchwenn et all ( 20051) found the respective 
rates of false and missing alarms for this correlation to be 15 and 20 percent, 
which underlines the need for more sophisticated approaches. A first correction 
is the tendency of the Parker spiral magnetic field to cause a slight westward 
deflec tion for fast CME s, and a stronger eastward deflection for the slower 
ones ( Wang et all 120041) . From a sample of 841 CME observations using the 
Large Angle and Spectrometric Coronagrap h (LASCO) on board the Solar and 
Heliospheric Observatory (SOHO) mission. 1st. Cvr et all ( 20001 ) found 14% of 
them to exhibit clear indication of non-radial motion. 

Obscrvationally derived trajectories are often ambiguous because coronagraph 
images suffer from projection effects, and the frequencies of the type II radio 
bursts associated with an eruption only provide the source heliocentric distance 
(which furthermore depends on the assumed electron density profile). Corona- 
graph data fro m the Solar Terrestrial Relations Observatory (STEREO) mission 
( Kaiserl | 2005[) can be used for stereoscopic reconstruction from image pairs 
( Howard and Tappinl. 2008). Thi s, however, is a well-defined problem only for 
curve-like objects (jlnhesterl . 12006 ). while for extended, diffuse objects like CMEs, 
some ambiguity remains. This is usually resolved by manual identification of 
those points in image pairs which are assumed to belong to the same point-like 
feature in physical space; since this procedure uses two 2-D points to determine 
a singe 3-D poi nt, it is in fact an over-defined pr o blem. For the small sample 
investigated by iMalonev. Gallagher, and McAteerl ( 20091) . these authors found 
their derived trajectories to be consistent with quasi-radial expansion. 



2.2.2. Simple kinematic models 

Several authors have tried to deduce a CME's travel time T in which it covers 
a distance R s / C as a function of its initial velocity vq by fitting formulas like 



v T + aT 2 /2 = R s/C 



(1) 



to (vq,T) data pairs from actual events, thus using the required acceleration 
a (due to thermal pres sure, magnetic forces, gravity, momentum conservation 
("snow plough effect", Tappin . 20061 ) . and aerodynamic drag) as a free pa- 
rameter. iGopalswamv et al\ ( 20011 ) found a minimum variation of AT s=s 10 h 
for their sample if a = beyond 0.7 5 AU. The tentative inclusion of a drag 
term a(v) oc (v — v sw ) 2 ( Cargilll . 12004 ) was shown to have little influence on the 
method performance, though it is indeed observed that for larg e distances, the 
CME speed will approach the speed v sw o f the ambient medium (JLindsav et al 



19991 : IMalonev. Gallagher, and McAteerl . l2009h . Relating these empirical T 
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T(vq) formulas to observed events yields large discrepancies, which can be at- 
tributed to the inherent oversimplification of this method. Therefore, it has to 
be concluded that both the CMEs themselves and the interplanetary medium 
which they encounter are much too variable and structured for simple fitting 
laws of this kind to yield more than rough estimates. 

As an intermediate step between one-parameter fitting and fully self-consistent 
MHD simulations, spatia l ly resolved kinem atic models such as the HAF model 



(JHakamada and Akasoful . Il982t iFrvl . 119851) have been used to predict the ar 



rival times of shocks and the spatial structure of the inner hcliosphcre as it is 
modulated by corotating interaction regions. This model projects a (possibly 
time-dependent) boundary condition at 2.5i?© outwards along strea m lines, and 



has b een calibrated using direct (albeit only 1-D) MHD simulations (pun et all 
19851) . 



2.2.3. Numerical MHD propagation models 

The numerical study of propagating CMEs is a challenging task mainly because 
the need to track structures with details smaller than a solar radius across at 
least 1 AU requires very high spatial resolution. On top of that, the configu- 
ration has no apparent symmetry properties that could be exploited to reduce 
computational expenses. During solar maximum, the Sun's global magnetic field 
is inherently 3-D, and even during solar minimum, CMEs arc never observed to 
travel along the Sun's polar axis, thus breaking the rotational symmetry of the 
field. 

In the past, two strategies have been used to deal with this difficulty. The first is 
to ignore the misalignment of expansion direction and magnetic sym metry axis 



and to assume cylindrical symmetry anyway {e.g. IChane et all 120061 ) . This im- 
plies a closed, torus-shaped CME geometry which is not anchored on the Sun, a 
configuration that could potentially be releva nt for tube-shaped magnetic clouds. 



In the ir comparison of 2-D versus 3-D models. iJacobs. van der Hoist, and Poedts 



(2007)) conclude that propagation models with cylindrical symmetry, however 
inadequate from a principal point of view, can still provide useful and computa- 
tionally inexpensive estimates, which can then be used to set up a refined and 
symmetry-free follow-up simulation. The general usefulness of such approaches 
requires that the parameters are properly transferred between both geometries, 
which necessarily requires several ambiguities to be resolved. 
The second avenue is to use a fully 3-D model with sufficient resolution in all 
three coordinate directions. The resulting vast increase in computational expense 
can be moderated by the use of specially tailored grids, in particular spherical 
grid geometri es with a radial mesh spacing Ar which increases with heliocentric 



radius r (e.g. iPomoell. Vainio. and Kissmannl 120111 ). 



For a global fluid simulation of CMEs, the numerical grid must be chosen 
carefully to optimize the trade-off between the advantages and shortcomings 
of different grid geometries. While spherical, Sun-centered grids are obviously 
best adapted to the Sun's shape and thus greatly simplify the specification of 
boundary conditions and the description of predominantly radial expansion, the 
large variation in cell sizes may lead to undesirably low time steps, and the 
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necessary inclusion of the polar axis requires a delicate treatment of coordinate 
singularities. The latt er can be avoided by the use of so-ca lled overset grids, such 
as the Cubed Sphere (IRonchi. Iacono. and Paoluccil 1 1996) or the Ying- Yang grid 
(JKagevama and Satol . 120041 ) . These, however, require some form of interpolation 
scheme to establish a seamless connection between the sub-grids, which usually 
destroys the conservative properties of the code within the int erpolation r egion . 
For a pplications of these grids to 3-D solar wind modeling see lFeng et all ( 20101 . 
20111 ). Cartesian coordinates, which do not present such difficulties, remain 
popular among CME modelers for exactly this reason, although the implemen- 
tation of spherical boundaries like the photo sphere then either req uires some 
form of weighted interpolation procedure (e.g. iKleimann et aZ.I . I2009T ) or a much 
increased spatial resolution. The latter requirement becomes m uch less severe by 



the u se of adaptive mesh refinemen t techniques (e.g. BATS-R-US.Ide Zeeuw et al.. . 
20001 ) . Logically rectangular grids (jCalhoun. Helzel. and LeVequel . 120081 ) , which 



permit the smooth inclusion of curved boundaries into an otherwise Cartesian 
grid geometry, harbor the potential to combine the advantages of Cartesian and 
spherical coordinate grids, but have until now not been used in any large-scale 
3-D simulation with relevance to space physics. 



Anot her interesting approach is the use of multi-scale models ( e.g. iRilev et al 



20061 ). which consist of a radially nested sequence of grids, such that the outer 
boundary condition of a given grid is fed into the next, coarser grid as an inner 
boundary condition of the latter. The fact that the solar wind plasma flow 
becomes super- Alfvenic after only a few Rq frees these authors from the require- 
ment to run the model simultaneously on all sub-grids. The multi-scale approach 
is also a key feat ure of the Space Weather Modeling Framework (SWMF, see 
Toth et all . 120051 ) . which condenses a total of nine model components (for the 
magnetosphere, inner hcliosphcrc, radiation belts, etc.) into a modular, versatile 
space physics-based simulation environment. 

2.3. Classification by aim and scope 

Existing models for the expansion/interaction phase of CMEs largely fall into 
two categories: 

1. 'principal' studies, which intend to investigate the relevance of a (usually 
small) set of free parameters on the resulting dynamics, and 

2. 'realistic' simulations, which are aimed at a usable prediction of the devel- 
opment of actual events, and therefore need to include as much physics and 
data as possible. 

Consequently, these two types of models differ noticeably in the way in which 

they implement boundary and initial conditions. On the one hand, principal 

models tend to use a backgroun d solar wind which is either uni form fe.q.lVandas. Odstrcil. and Wataril | 
2002 c iDalakishvili et all . 2011b ) or moderately structured (e.o. lOdstrcil and Pizzol .[ 
1999c iManchester et aUl2004bl ). combined with a dipolar, or at most quadrupo- 
lar magnetic field, and possibly amended with a planar current sheet. Realis- 
tic global-scale models, on the other hand, routinely rely on synoptic magne- 
tograms to specify a potential coronal magnetic field as a photosphcric bound- 
ary condition, from which a reasonable initial condition can be derived (e.g. 
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Havashi. Zhao, and Liul . l2006t iLugaz et all 12007 ). To fully constrain this field 



solution, additional assumptions about the type of B field need to be made, the 
most popular being that the field is either potential (V x B = 0), linearly force- 
free (V x B = qB), or non-linearly force-free (V x B = a( r) B). For a review o f 
strategies and algorithms to obtain force- free solutions see lMetcalf et all < 20081 ); 
a critical asses s ment of their performance on actual data has been presented by 
iDe Rosa efal\ (J2009f) . 

Alternatively, vector magnctograms have been used to calculate localized non- 
linear force-free fields on the scale of active regions, f r om w hich CME-like erup- 



tions can be launched in a simulation ([Kataoka et all . 120091) . Once the magnetic 



field is known, empirical formulas can be used to determine suitable initial con- 
ditions for the remai ning fluid quantities density, velocity, and temperature (e.g. 
Detman et oil . l2006h . 



2.4. Method of initiation 



Given that the physical process which triggers an eruption is still not conclu 
sively settled, there is some freedom in the choice of methods to numerical! 
start this proces s . Den sity-driven models (jGroth et all l2000t IChane et all 1200' 



Kleimann et all . 12003 ) increase the pressure and/or density (and occasionally 
also the radial momentum, th ereby imposing an additional mass inflow (e.g. 
Keppens and Goedbloedl . 120001 )) below a closed field line configuration until the 



inward tension of the field can no longer confine the growing internal pressure 
and gives way to a violent expulsion of the accumulated plasma blob. Similarly, 
addit ional magnetic f l ux can be artific ially pushed through the photosphere 



(e.g. iFan and Gibsonl . 120071 : iFanl . 12011 ). Note that this latter approach must 



be carefully distinguished from flux emergence models, which arc used to study 
the physical process of buoyant flux tubes emerging from the convection zone. 
The initial field may be furthe r destabilized by, e.g . , loca l magnetic reconnection 
(e.g. iForbes and Priestl . Il995t IChen and Shibatal . 120001 ) . or by an initial force 
imbalance, e.g. due to buoyancy. Popu lar initial setups o f the la tter kind include 
the self-similar ex panding flux rope b y Gibson and Low! (J1998I) , which has been 



employed by, e.g.. Manchester et all (J2004bl ) t o initialize simulation s of C MEs, 



just like the famous Titov-Dcmoulin flux rope ( Titov and Demoulinll 19991 ). 
Finally, the required instability can also be produced by shearing the photo- 
spheric footpoints of an arcade until it becomes unstable and is forced to open up 
spaccwards. The effect of s hearing mot ions on solar mag netic structu r es has been 
investigated the oretically (ILowl. Il977l). observationallv (IDeng et all. 120011). and 



numerically (e.q. lMikic and Linken . il 994 Uacobs. Poedts. and van der Holstl . 120061) .[ 



From these methods, the density-driven one is special in that it is not supposed 
to mimic a process that is actually assumed to happen on the Sun; moreover, its 
justification lies in the pragmatic notion that if the focus is on the propagation 
phase, the details of initiation are of secondary importance, thus justifying the 
use of a method that is easy to implement an d produces a relatively realistic 
outcom e. The same rationale was also given bv lPomoell. Vainio. and Kissmann 
(|2008l ). who start their Cartes ian 2-D simulation using a set of line currents 
similar to the one employed by IChen et all ( 20021 ) , and then apply an artificial 
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volume force to the interior of the detached flux rope to pull the core region up- 
wards (which requires the trajectory of interior fluid elements to be dynamically 
traced during the acceleration phase) . 

2.5. Treatment of the energy budget 

Despite almost six decades of undiminished efforts to understand t he heating 



of the solar corona , this vital issue is still far from being settled (jCranmer 



20021 ; iMarschl . 120061) . For practical reasons, models seeking to describe a CME's 



propagation in the corona and inner heliosphere have to adopt a formalism to 
model the distribution and release of energy, which opens up another possibility 
to classify them. 

The simplest one would be to link the plasma pressure p and mass density p via 
an adiabatic closure relation 

p ~ f . (2) 

In this case, again the simplest (albeit also the most unrealistic) choice would 
be to assume an (almost) isothermal plasma by setting the adiabatic index 7 to 
unity, or a value slightly larger than unity, such as 1.05. Using the physically 
appropriate value of 5/3 in conjunction with Equation ([5]) is usually not feasible 
because the effect of adiabatic cooling would quickly reduce the temperature 
down to unrcalistically low values. This problem is sometimes mediated by using 
a spatially varying a diabatic exponent, either as an explicitly prescr ibed depen- 
dence 7 = 7(r) (e.g. iFahr. Bird, and Ripkenl . Il977t IWu et al. ._ 1993). or derived 



from energy considerations based on the Bernoulli equation (jCohen et all . 120071 ) 
or from interpreting the low value of 7 close to the Sun in terms of internal 
energy being stored as additional degrees of freedom in large-scale turbulence 



( Roussev et all . 120031 ) . The obvious downside of using unphysically low adiabatic 



indices is the resulting inability of the model to correctly reproduce the properties 
of shocks, which may or may not be acceptable in a given situation, but which 
must in any case be kept in mind when interpreting the results thus obtained. 
If 7 = 5/3 holds throughout the considered volume of space, a full energy 
equation for the total energy density e like 

a t e + V-[(e+p+||B|| 2 /2)u-(u.B)B] = (g • u + Q)p (3) 

(in which u denotes the flow velocity and g the acceleration due to gravity) with 
heating source term Q has to be used. Popular choices for Q include ad hoc 
heating functions like 

Q(r) = q(v) [T - T] (4) 



([Manchester et all l2004b| ). where To is a spatially dependent "target temper- 
ature", which is fine-tuned to reproduce a realistic quiet-Sun temperature dis- 
tribution. The main shortcoming of this approach is that it is obviously biased 
towards the target temperature, or in other words, it will systematically underes- 
timate the temperature change which is brought about by, e.g., the passage of a 
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CME. Another, more recent approach is that of lPomoell. Vainio. and Kissmann 



(|201ll ). who employ an energy equation 

d t (p/pi) + vi-V(p/pi) = S (5) 

and first run a quiet-Sun simulation using 7 = 1.05 and S = (and no CME), 
followed by the 'main' run using 7 = 5/3 and 

S = ui-v(jn/p$ /3 ) (6) 

where the subscript 1 denotes values from the 7 = 1.05 run. This procedure 
guarantees that both runs share the same stationary (quiet-Sun) state. 

2.6. Alfvenic wave heating 

As an alternative to the pragmatic heating functions of Section 12.51 attempts 
have been made to amend the underlying solar wind model with a more physical 
heating process. A promising candidate for such a mechanism is heating by 
Alfven waves, which are generated at photosphcric levels, then travel outward 
along magnetic field lines while being shifted towards an upper limit frequency 
/h at which they dissipate and deposit their energy as heat. In its full form, this 
scheme requires an entire wave power spectrum P(f, r, t) to be modeled for each 
position r. A dynamic equation for the temporal change of P has to be solved for 
each frequency /, and the waves then couple back to the MHD system via: 

1 . the accelerating negative gradient of the wave pressure 

Pw ~ I h P(f) d/ (7) 

Jfo 

representing the magnetic pressure of the fluctuations, and 

2. a wave heating term 

Q w = F(f h , r) - P(/ h (r),r) [u ± v A ] • V/ h (r) (8) 

in which the so-called cascading function F = F(P, /) governs the spectral 
evolution of P, and va is the local Alfven velocity. 

More details can be found in the review bv lFichtner et all ( 2008 ). This formalism 
has been successfully applied to purely radial models of the quiet-Sun solar wind 
( Laitinen. Fichtner. and Vainiol 120031 ). and appears to be a promising alternative 



to the phenomenological heating functions discussed in Section 12.51 although 
the extension to the generally more complicated magnetic field structures en- 
countered in symmetry-free CME simulations is clearly non-trivial. First results 
obtained with a wave-heated solar wind model (albeit relying on the scalar wave 
energy densities e± parallel (+) and anti-par allel (— ) to B rather than the full 
spectral information) have been presented bv lvan der Hoist et al\ ( 20101 ). 
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3. Connecting to observations 

As can be expected from their very different aims, 'principal' and 'realistic' 
models give qualitatively different types of results; one could also say that they 
provide answers to different types of questions. 

3.f . Selected findings from 'principal' models 

Since several analytical CME models rely on the simplifying assumpt i on of 
self-similarity (IFarrugia. Osherovich. and Burlagalll995tlGibson and Lowl .lf998; 



Nakwacki et qZ.I , l2008tlWang. Zhang, and Shem 20091). it is a vita l questi on whether] 
this property can be confirmed numerically. iDalakishvili et all ( 2011al) consid- 
ered an idealized model for a cylindrical magnetic cloud, and numerically con- 
firmed that an initially self-similar configuration maintai ns this property a t later 
stages. Using a slightly more realistic 3-D CME model. iKleimann et al\ ( 20091 ) 
compared the derived time profiles of density and magnetic field strength at dif- 
ferent heliospheric distances, and also found indication of self-similar evolution at 
least in the early phase of propagation, which in this case covered a mere f Rq. 
This is consistent with obser vations which show a CME's cone angle to be ap- 
proximately constant in time (jSchwenn et all 120051 ) , although other simulations 
have been performed in whic h the evolution starts to depart from sel f -simil arity 

(e.g. iManchester et all l2004at IChane et all I2006J) . As 



after several tens of R 



was demonstrated bv lRilev and Crookerl (J2004I ). these findings can be reconciled 



by observing that the passive advection of a structure in a spherically diverging 
flow causes the aspect ratio (azimuthal over radial extent) of the structure to 
increase linearly with radial distance. A CME of initially spherical cross-section 
will thus tend to be flattened into a pancake-like shell, a process counteracted 
only by the CME's own pressure-driven expansion and its magnetic tension force, 
which acts towards a reduction of field line curvature and thus tends to work 
against exces s ive fla ttening . 

Jacobs et all (J2005I ) compared several popular CME propagation models and 
found the resulting dynamics to be strongly dependent on both the included 
physics (most notably the heating mechanism) and the background wind, among 
other things confi rming the in t uitive notion that higher speeds are found in 
fast, dilute winds. IChane et all (|2006l ) investigated the influence of the CME's 
initial polarity with respect to the background field and found that 'normal' 
prominences result in fast, compact, approximately circular CMEs which get 
deflected equator wards, while CMEs which develop from 'inverse-' prominences 
are slower, pancake-shaped, and tend to deflect polewa rds, thus confirming a 
previous theoretical prediction by I Zhang and Low! ( 2004 ). 



3.2. The situation for 'realistic' models 

A functional prediction tool as envisaged by the space weather forecasting com- 
munity would take the observed distribution of photosphcric magnetic field 
at the instant of CME outbreak, infer from it the global initial state for all 
MHD variable fields, then run a detailed CME simulation from which the ex- 
pected physical conditions at a given location (typically at Earth orbit) can 
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be extracted. Therefore, realistic models will typically try to reproduce in situ 
observations by spacecraft such as WIND, Ulysses, or the Advanced Composition 
Explorer (ACE) (see, e.q.. lChane. Poedts. and van der Holstl . 120081 ). For regions 



closer to the Sun, another possibility to assess the predictive capabilities of the 
model is the creation of artificial white-light coronagraph images, which can 
be contrasted with actual images of the respective event, as was done by, e.g., 
iLugaz et all ( 20071). A first quantitative comparison of this type was presented 



by iManchester et al\ ( 20081) . who not only compared the absolute projected 



brightness distribution but also the derived mass and velocity of the CME in 
question. 

While visual inspection of these comparisons often indicates impressive agree- 
ment between simulation and observation, it must however be noted that such 
comparisons tend to be biased towards good agreement, since publications will 
typically show only a few selected results from a limited parameter range, in 
spite of the fact that models are quit e sensitive to the chosen paramete r s (e.g . 
ISchriiver et all . 120081) . (For instance, IChane. Poedts. and van der Holstl ([20081 ) 



perform 24 different runs but explicitly state that they only show their best 
fit to the ACE data of the target event.) This makes it difficult to assess the 
predictive capabilities of a given model for 'real world' forecasting applications, 
in which a posteriori selection of parameter sets is not a valid option. 



4. Conclusions 

CMEs represent a very diverse and important class of heliospheric transients. 
Besides a purely scientific desire to understand their true nature, the urgent 
need for accurate space weather forecasts creates a strong commercially-driven 
incentive to intensify the study of these phenomena. 

Since purely kinematic models are inadequate to capture their highly non-linear 
evolution, 3-D MHD simulations are indispensable to model the different stages 
of a CME's life cycle, with the long-term goal of providing reliable forecasts. 
Associated difficulties mainly concern the need for high spatial resolution - 
although the ever-increasing performance of computing power has helped much 
to mitigate this problem — , the ambiguity induced by the need to complement 
observational input by reasonable assumptions about the remaining quantities 
not obtainable from direct observation, and the fact that to this day, no univer- 
sally accepted model for the initiation phase exists. 

Generally speaking, the model developing community draws huge benefits from 
high-quality spacecraft data to: 1. constrain the initial and boundary conditions 
and 2. to allow for an a posteriori validation of the obtained simulation results. 
Since these results and predictions crucially depend on both the adopted initial 
parameters and the physical effects included in the model, more (parameter) 
studies are needed to quantify their respective influences. Furthermore, an un- 
biased comparison between different space-weather prediction tools is currently 
difficult because each group of authors tends to choose a different event, relies 
on different input data from the event in question, and sometimes also picks 
different methods and criteria to validate their results. It thus seems that the 
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community could benefit enormously from some form of standardized space 
weather forecasting benchmark, in analogy to community-wide b enchmark test 



suite s known from and used for the simulation of the heliq s phere (jMuller et al. 



2008 ) or for convection in the Earth's interior ( King et aZ.I . I2010TI 
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